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Abstract 

We are presenting an algorithm capable of simplifying tensor polyno- 
mials with indices when the building tensors have index symmetry prop- 
erties. These properties include simple symmetry, cyclicity and those 
due to the presence of covariant derivatives. The algorithm is part of a 
Mathematica package calle d Tools of Tensor Calculus (TTC)[web address: 
[rttp://baldufa.upc,es/ttc| 

1 Introduction 

The two main languages commonly used in writing tensor calculus expressions 
are the intrinsic notation and the index notation. The intrinsic notation seems to 
be the preferred language for making computer algorithms work with symbolic 
tensor expressions. The index notation, on the other hand, is extremely powerful 
in expressing and manipulating tensors. There are some simple expressions 
(in index notation) which are difficult to express in intrinsic notation for the 
purpose of introducing new operators with new properties. This is the case e.g. 
expressing the simple operation of raising an index ([[J) 

V k =9 3m T lmk (1) 

Recently some computer tools capable of working with index notation at a 
symbolic level have appeared 0, D • 

In the present study we are facing the problem of tensor polynomial simpli- 
fication using index notation. In a recent article || we dealt with this problem, 
restricted to the case of simple index symmetry properties of the building tensors 
like 

T Ml) ... M „ ) +aT ll ... 4 „=0 (2) 

a being a scalar and a a permutation of {1, n}. (||) includes symmetry and 
antisymmetry of any sequence of indices, and pairsymmetry. We will refer to 
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these kinds of properties as monoterm. The algorithm proposed in [|[ does not 
use any kind of library, either temporal or permanent, and for that reason it 
does not require a lot of computer memory. Nevertheless, beside the monoterm 
restriction the algorithm presented did not draw benefit from all the calculus 
done. It was reasonably fast on one monomial but not on a hundred, which is 
a problem to be solved when using other properties. 

V.A. Ilyn and A.P.Kryukov have proposed a method capable of working 
with index symmetry properties, which include cyclic type properties like 

Tijk + Tkij + Tjki = (3) 

was included by the authors in a set called multiterm like properties. Al- 
though the problem was theoretically satisfactorily solved, the computer time 
and the amount of computer memory needed to make elementary simplifica- 
tions was large enough to consider that this study may have some practical 
and fundamental limitations ( since an expression with 11 indices needs about 
2500Mb). 

One of the problems encountered by V.A. Ilyn and A.P.Kryukov is related to 
the dummy indices. They consider the possibility of renaming dummy indices 
as new relations to add to the set of relations generated using the true index 
properties of building tensors. So in the monomial Sj^T^ 1 as well as the prop- 
erties due to the fact that S is a symmetric tensor they add the ones deduced 
from the fact that the indices j, k are dummy, i.e., 

S jk T kjl - S k] T 3ki = (4) 

The number of this kind of relations increases factorially with the number of 
dummy indices, so the algorithm needs to work very hard for monomials of 
interest. Other authors ||, || have studied simplification algorithms but these 
need to be implemented. R.Portugal Q handles cyclic like properties similar to 
|3) . We think that his algorithm can work very well for monoterm properties ( 
like H does) but not for multiterm ones. The trouble is that although apparently 
he can decide if a multiterm rule takes advantage using them increasing the 
lexicographic order of the working monomial, it can happen that this decision 
must be delayed after analyzing the effect of some chain of multiterm rules 
applied to the monomial. In general the whole set of monoterm plus multiterm 
monomial equations must be considered. 

The algorithm that we present here, written in Mathematica Jj6| language, 
is a part of the package Tools of Tensor Calculus (TTC) [Q, ||, |^, from now 
on we will call both, the algorithm and the package, TTC). TTC can work 
with monomials whose building tensors have both monoterm and multiterm 
properties. In addition to dealing with dummy indices it also handles symme- 
try (including the partial derivatives, although it is not recommended because 
of their non-tensor character and the increase of the corresponding computer 
time. In this paper we will leave aside partial derivatives ) and antisymmetry 
of any sequence of indices, pairsymmetry, cyclic type properties and the Ricci 
commuting relations due to the presence of two or more covariant derivatives. 
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Taking the Riemann tensor as a paradigm TTC can handle monomials that 
have this tensor as a factor, among others, with the known properties: 



Rijkl = —Rijlk 


(5) 


Rijkl — Rjikl 


(6) 


Rijkl = Rklij 


(7) 


Rijkl + Riljk + Riklj = 


(8) 


Rijkl-.m ~t~ Rijmk\l ~t~ Rijlm;k 


(9) 


.. ;m + (Riemann terms) .... = 


(10) 



fe-2 



where ( |Io|) refers to Ricci commuting relations. We have left the inclusion of 
the properties derived from the dimension of the base space, for the Riemann 
tensor and for the generic case, for a future work [ fio| . 

We have applied the algorithm presented in a real calculus ( superenergy 
tensor problems) finding new relevant results Jill. 

The plan of the present paper is as follows: in section 2 we present the 
problem and explain our solution from the theoretical point of view. In section 

3 we describe how the algorithm deals with Ricci commuting relations. Section 

4 is devoted to describe the main function of the algorithm implementation. In 
the Appendix A we show some examples of how TTC works in practice. 



2 The simplification problem 

The problem can be stated as follows: 

Given a tensor monomial written in index notation together with all index 
properties of the building tensors, to find if it is zero and, if not, a unique 
canonical equivalent polynomial. 

The canonical equivalent polynomial is derived from the complete set of 
monomial relations generated from the exhaustive application of all index prop- 
erties of the tensor factors over the original monomial. The solving procedure is 
related to some ordering criterion of the set of monomials. The concrete order- 
ing criterion is not important, we will assume that we have a unique ordering 
criterion. 

Free and dummy indices 

Let M(Ip, Jf)) be a given tensor monomial, with F free indices Ip = {jpm, ■■■,jp(F)} 
and D dummy indices Jf) = {j q (i)j ■■■>j a (D)} j ct an d (3 being any permutation 
of {1, ...,£>} , {1,...,F} respectively. Example: 

M({i,k,l},{j,m,p}) = S llk Ti p m T p lrn 
Due to the meaning of the dummy indices the following identities hold 
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M(4,J«) = M(4,J D ) Va (11) 

Jd being the canonical ordering of dummy indices in the monomial. This means 
that there will be many equivalent expressions for a given monomial. This is 
the first problem to be solved. 

In H we looked at a method in which a monomial could be rewritten in a 
canonical form with respect to the dummy indices independently of the num- 
ber of dummy indices present in the monomial, without using much computer 
memory. 

Essentially the method consists of using an ordering criterion that can be 
applied to a monomial without prior consideration of all the monomials to be 
ordered. The method presented here is an improved version of ||, in order to 
optimize the amount of storing data and computer time, following the criteria: 

1) Free indices are renamed: IndexF[l] , IndexF[2] , . . . all of them in 
upper position ( contravariant) 

2) Dummy indices are renamed IndexD[l] , IndexD[2] , . . . all of them 
in upper position ( contravariant) 

3) We find the permutations of free, f3, and dummy, a, indices in order to 
find the most ordered monomial expression. Note that at this stage the result is 
independent of the original order of free indices. We call this globalcodification: 

M(4, J a D ) ^ M lcatl ° n M(I F , J D ) = M(I F ) (12) 

being Ip the canonical ordering of free indices in the monomial. 

M(If) is the skeleton of the monomial. M(If) carries all the necessary and 
sufficient information on the monomial independent of the name of dummy and 
free indices. Different monomial skeletons will be denoted by Mi(I F )M2(I F )... 

Finally we take the permutation f3 performed on free indices over the global- 
codification of the monomial. The final result is called codification of the original 
monomial and the process will we called (O)-rules. 

M(4, JS) (0) ^' CS M (4, J D ) = Jf (»>(/£) (13) 

Note that the canonical sets If and Jd are defined through ( [l2|) since they 
depend on the skeleton of the monomial. Also note that any permutation of 
free indices of a globalcodification corresponds to a codification version of the 
corresponding permutation of free indices in the original monomial, so when a 
monomial has many free indices this procedure optimizes computer time and 
memory. 

On the other hand, as we will see later, all properties found during a cal- 
culation are stored in codification version but all monomials with the same 
globalcodification version can use, due to what we have stated on the precedent 
paragraph, the corresponding /3-like permutation version of the stored proper- 
ties, so the work of finding properties is essentially independent of the order of 
free indices. 

In the following sections we will handle monomials in the form M^\4)- hi 
this way we will not mention the dummy indices problem anymore. 
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Tensor properties 

Of all tensor properties we will emphasize those that will be of importance 
within the algorithm. 

• Monoterm properties: the result of their application over a given mono- 
mial is also a monomial. 

• Multiterm properties: the result of their application over a given mono- 
mial is a polynomial (more than one term). 

• Homogeneous properties: they are properties of only one tensor, say 
^lii... i n 3 °f t ne form 

M 

T ll ..,„+^a CT3 T Vj(1) ...^ ( „ ) =0 (14) 
i=i 

<jj being permutations of indices and a„. scalars, in such a way that ( pTj ) 
is invariant under any of the ctjT 1 , that is 

M 

j=l 

M 
3=1 

^4^ being some scalars. Example: see (J3j) . 

• Inhomogeneous properties : These are properties relating one or more 
tensors whose form is not reducible to (|l~4| ) fulfilling ([l5|) . 

The most important example of inhomogeneous properties is the Ricci com- 
muting relation. The simplest case is when we have a vector field v, then 

Vi;j;k ~ Vi-k;j ~ V m R^j k =0 (16) 

R being the Riemann tensor 

The advantage of handling homogeneous properties is that we can pick up 
a term and use it as a pattern to be applied when finding rules over a given 
monomial. This procedure, due to the homogeneous character, will be equiv- 
alent to use all terms in substitutions on the monomial. This is not true for 
inhomogeneous properties and special treatment must be considered. 

On the other hand some inhomogeneous properties can be applied directly 
iff the right hand side inherits all the properties carried from the left hand side. 
There is an important example when defining the Ricci tensor. R"l m j = Rij can 
be applied directly iff Ricci inherit from Riemann all the resulting properties, 
for example 

R%n = \RZ, (17) 
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which is a consequence of applying Ricci definition to the original Bianchi II 
relation. If we want to use the definition of the Ricci tensor as a direct relation 
we must include also as a direct relation, as well as the deduced properties 
belonging to the Ricci tensor itself, Rij—Rji = 0. Although such a procedure can 
be difficult and unnecessary because we can work exclusively with the Riemann 
tensor, it is very convenient in order to eliminate superfluous indices. 

Induced monomial properties 

Given a set of monomials M = {M± (I F ), M 2 \lp), ...} and a set of properties 
P = {pi, ...pk--} of the building tensors, the set of monomial properties induced 
by P, P(M), is the set of properties which arise applying P over M and over 
all new monomials generated until no new monomials arise. 

The essential idea when simplifying a given monomial M^(I F ) is to take 
the set of properties of the building tensors, P, and generate the set of monomial 
properties P({M^(I F )}), solving them with respect to the ordering criterion 
and finally applying the corresponding rules onto the original monomial. The 
trouble is that the set P({Af(°> (/£)}) can be too large to be reasonably gener- 
ated and solved. 

It is interesting to try to break the set of properties P into two sets, say 
Pi and P 2 . Then generate and solve P\({M^ (I F )}) applying the resulting 
rules over M^(I F ) giving, in general, a polynomial Pol^ built by the set of 

monomials {..., M^\l^), ...}. Then generate and solve P 2 ({..., M^ 1 \l F ), ...}) 
to give a new set of rules to be applied onto the polynomial Polm giving the 
final result. 

We will call separable properties the set P\ of P which can be separated 
from the rest in the sense of the above paragraph. The separability is a possible 
property of the given set P which must be demonstrated case by case. Properties 
which affect different tensors or different sets of indices in the same tensor are 
clearly separable, but this is not only possibilitu. It is easy to assess whether a 
given monoterm property is separable or not. The most important case is the 
Riemann tensor 

Lemma 1 From the whole set of properties of the Riemann tensor, Rijki + 
Rijik — and Rijki + Rjiki — are separable. 

Proof: We use the Riemann properties P = {(g), (g), @, (§), (§), @}. Let 
Pi = {piM = {§), (§} and P 2 = P - P x . 

Given any monomial having the Riemann tensor as a factor and using Pi to 
improve the ordering, can be written as 

RF(i,j,k,l) = ...Rijki... 

where the indices ijkl can be free or dummy but we can consider that the 
actual names will correspond to the most orderly form. Pi are separable iff 
P 2 ({RF(i,j, k, I)}) = P 2 ({RF(i,j, I, k)}) = P 2 ({RF(j, i, k, I)}) 
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Let us analyze P2 case by case. We will use the notation (re/.)({M}) to 
indicate the set of properties Pi (M) being Pj a property of the set P: 
pairsymmetrie (|7|) 

Clearly {(g) M} is separable from the set {@, (|), (|)}. 
Bianchi I (p) 

{§({RF(i,j,k,l)}) gives 

RF(i,j,k,l)+RF(i,l,j,k) + RF(i,k,l,j) = (18) 

(i ({-R^CU, fc)}) gives RF(i, j, I, k)+RF(i, k, j, l)+RF(i, I, k, j) = which 
after using 

{(|),(§)} coincides with (|T|) 
(§({RF(j, i, k, OV) gives RF(j, 1, k, I) + RF(j, k, I, i) + RF(j, I, i, k) = which 
after using {(||), @} gives 

- EF(i, j, k, I) - RF(j, k, 1, 1) + RF(j, I, i, k) = (19) 

which does not coincide with (|lg| ) so ^2 is not separable from Bianchi I, but since 
we have also pairsymmetrie we can consider (||)({i?F(/c, I, i, j)})which gives ( [l9| ) 
so we can conclude that {(§), f§} is separable from the set {(|), (|), (0), (|)}. 
Bianchi II @ 

Considering a monomial of the form RF(i,j, fc, Z, m) — Rijki-m---- and finding 

©({-ft-PXij h I, m)}) and ^)({RF(i, j,l,k,m)}) the properties generated are 
the same as the property generated by M)({RF(i, j, k, I, m)}), so {(||), (||)} is 
separable from the set { (|) , (|) , (|) , (§) , (|) } . 
Ricci commuting relations dlOj ) 

It is easy to observe that the (Riemann terms) of ( |l0| ) have the same symme- 
try properties with respect to the indices i, j, k, I as the term Rijki-m-m specially 
{(§, (§}, so clearly {(f), (§} is separable from the set {(g), (g), (§), (|, ®, ©}, 
and this completes the proof. 

T"TC uses this separability property when internally defining the Riemann 
tensor. 

In general we assume two sets of properties accordingly to the possibility 
of separating some from the whole set. We will call them (n)-properties with 
n = 1,2 

How to handle (ri)-properties 

Here we explain how TTC finds, solves and applies (l)-propcrties and (2)- 
properties. We will supose that we depart from a monomial M^^Ip 1 ) so let 
n = 1 and go to step n. 

Step n: 

Given a monomial M[ n ^(Tp 1 ) the exhaustive application of the (rt)-properties 
generates a set of relations of the type 
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Ml n - 1) (4")+J2b ijpq M^- 1 \l^)=0 i,j = l,... (20) 
jq 

where i,j,p,q enumerate possible M, /?, b objects appearing in (^p| ) and b... are 
scalars. 

This set of relations can always be solved in favor of the subset of the most 
orderly monomials 

Mt- 1 \4n = J2 c ^ A t ) ( I F q ) (2i) 

j,q 

c ipjq being scalars and M^ n \lp q ) the most orderly monomials from the set 
M ^-i)( 7 A-) of @ Thc specific 

way to solve (^0|) will be explained in section 

4. 

Finally we apply the obtained (n)-rules over the original monomial 

^ s J2cii jq M^(4") (22) 
j,q 

(n)-rules are conveniently stored for further use, as we will see later 
If ri < 2 set n = n + 1 and go to step n 

3 Ricci commuting relations 

As we have mentioned, Ricci commuting relations are inhomogeneous properties 
which must be handled with special care. If we want the algorithm to work 
correctly it is not sufficient to take into account Ricci commuting relation in the 
form (where we are taking an example of a single vector held v) 

Vi-j-k -> Vi-k-j + v m K m i3 k (23) 

since these rules will not benefit from other rules one might have in the 

case of monomials which do not have covariant derivatives explicitly, as in 

To solve this problem we have built a function, called AntiRicci which takes 
any monomial and gives all the monomials related to it by Ricci commuting re- 
lations via covariant derivatives. After that we apply Ricci commuting relations 
over the covariant derivatives. 

AntiRicci [v m R m l3 fe ] = {v,.^} (24) 

With respect to separability we can write 
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Lemma 2 Given a tensor Tijki... all symmetry properties of indices (i,j,k,l..) 
are separable from the Ricci commuting relations 

proof: Is trivial due to the fact that Ricci commuting relations are invariant 
under these symmetries. 

This Lemma does not apply for properties involving covariant derivative 
indices. There is an important example when defining Vi as Killing vector field: 
Vi-j + Vj-i — is not separable from Ricci commuting relations. 

The examples on the Appendix A are specially addressed these questions. 

Ricci commuting relations are internally included in the set of (2)-properties. 
Accordingly, if a tensor has no-separable properties and we want to use Ricci 
commuting relations properly, the set of properties must be included as (2)- 
properties 

4 Algorithm description 

The algorithm for simplification of polynomial tensors in index notation, called 
Simplif yAlllndex, has the following functions: 

Functions related to single tensors 

• InputTensor: 

InputTensor ITsymbol, basis, ranklisf] declares that the symbol Tsymbol 
will be a tensor in the basis basis and with rank and type ranklist. 

Example 

InputTensor [T , XX , {1 , 1 , 1 , 1}] 

• InputSymmetries \_Tsymbol\_index\ ,symmetriespecifications'] declares that 
the tensor, previously introduced using InputTensor, Tsymbol has the 
symmetry properties of the indices, positioned as indicated by the argu- 
ment index, according to the specifications symmetriespecifications. 

Example 1: 

InputSymmetries [T [i , j ,k,l] , {i ,k, 1} [1] ] declares that T is totally 
symmetric with respect to the indices {i,k,l} and add this property 
to the set of (l)-properties 

Example 2: 

here we declare the same properties as the Riemann tensor, see ||, 
||, [nj), for the tensor T accordingly to the separability theorem given for 
the Riemann properties. 

InputSymmetries [T[i, j ,k,l] , 

{{i,j»[l] ,{{k,l»[l] ,{{i,j>,{k,l»[2] , Cyclic [j,k,l] [2]] ; 
InputSymmetries [T[i , j ,k,l , . ;m] .Cyclic [k, 1 ,m] [2] ] 
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• BasicRules [n] : once the symmetries of indices have been declared the 
corresponding rules that implement these symmetries are stored in the list 
BasicRules [n] . n = 1 for (l)-rules and n = 2 for (2)-rules. The value 
n = is reserved for the case of properties that can be used directly , 
like T\ = 0, which are generally introduced explicitly by the user, or some 
properties TTC already has for Riemann tensor. 

Functions related to monomials 

In the following we will assume that there are no direct properties stored in 
BasicRules [0] as they are inessential. 

The main function is Simplif yAlllndex: 

Index [metricname,Slwpllf yAlllndex [2] ] [polytens] simplify, i.e. canon- 
ize, the polynomial polytens taking into account all the declared properties of 
the building tensor, as well as the dummy indices and the properties related 
to covariant derivatives. Simplif yAlllndex [2] is applied over each term of 
polytens. This function has, in turn, the following internal functions (We will 
assume that the initial monomial is M(I F , Jg)): 

• Simplif yAlllndex [0] : first we canonize the monomial with respect to 
the dummy indices. This step is taken each time a monomial is the result 
of some kind of index manipulation. So our initial monomial and all other 
generated monomials will have the form (I F P ) 

• LlistaMonomis \_label,n,counter(n) ] : in these lists all the monomials con- 
sidered before the present case are stored in globacodification version to- 
gether with its /^-permutation which relate it to the way they are present 
in the corresponding rules, n = 1 for (l)-rules and n — 2 for (2)-rules. 
counter (n) is a counter which says that the monomials therein are the 
counter (n)-th studied monomials. 

• LlistaRules \_label,n,counter(nJ] : in these lists the simplifying rules of 
the monomial in the corresponding LlistaMonomis are stored. 

The specific way TTC stores (n)-rules is 
LlistaMonomis \label,n,counter(nJ\ = j..., jM^™ -1 ^ (/p), , ...j 
LlistaRules Uabel,n, counter (nj] = j..., {m}™ _1) \f F v ) -> J2( 3 , q ) Koq^^^f )} » 

• Simplif yAlllndex [{n}] , with n = 1,2, simplify a monomial already 
simplified with respect to n — 1 properties, say M(™ -1 )(/p), with respect 
to (n)-propcrties. It searches if its globalcodification version, 

is in one of the corresponding lists LlistaMonomis and if it is the case it 
takes, if necessary, the corresponding (7 _1 o/3 fc ) permutation version of the 
LlistaRules solving with respect to the ordering criterion and applying 
the new rules to the monomial M^" -1 ) (Tp). 
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If M( n -^(I F ) is not found in any lists LlistaMonomis it applies the rules 
stored in BasicRules [n] exhaustively generating new lists LlistaMonomis 
and LlistaRules and applying them as explained above. 

• PFind is a simple tool that solves any system of linear equations with a 
fixed ordering criterion. 

Let {xi,X2, ...,xm} be the space of ordered variables to be solved. Let 
eqlist = {c 1 + a lj x 3 = 0, c 2 + ^ a 2j x 3 = 0, .., c N + ^ a Nj x J = 0} 

3 3 3 

(25) 

be the list of N equations with M variables ordered first taking into ac- 
count the variables and then the coefficients. That is to say, if in the u-th 
equation the most disordered variable is x m then x m +i cannot be present 
in the (n — l)-th equation. 

PFind takes the initial list of equations and starts the following loop: 

Step 1: PFind takes the first equation in eqlist. It is solved with re- 
spect to the most disordered variable. This result is stored in eqrule [1] . 
eqrule [1] is applied over eqlist and after eliminating zeros this result 
is assigned to eqlist. Go to step 2 

Step n > 1: If eqlist is empty the loop is finished. On the contrary 
PFind takes the first equation in eqlist. It is solved with respect to the 
most disordered variable. This result is stored in eqrule [n] . eqrule [n] 
is applied over eqrule [z] i = l,...,n — 1. eqrule [ri] is applied over 
eqlist, the zeros are eliminated and the result assigned to eqlist. Go 
to the step n + 1 . 

eqlist will remain empty if the system is compatible. These must be the 
case for a set of compatible tensor properties. 

The computer time and memory of PFind have been controlled. 

5 Concluding Remarks 

We have solved the problem of simplification of tensor polynomials in index 
notation with respect to a large set of properties for the building tensor satis- 
factorily. The computer time and memory needed are sufficiently reasonable so 
we think that the program can be of interest for practical purposes. We have 
successfully used the algorithms presented jll| 

Due to the generic character of the algorithm (the only requirement is that 
the input polynomial follows the usual mathematical rules) it can be useful in 
the kind of fields where tensors are a natural language. As is well known tensor 
calculus is powerful enough to cover large fields of mathematics, physics and 
engineering. 

There is an important set of properties not included in the present algorithm. 
These are the ones that appear when the dimension of the space is taken into 
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account. The essential reason for not including these relations is that applying 
properties of the building tensor do not generate this kind of relations. They 
arise when the whole monomial is taken into account so we need some tool to 
generate the dimensional properties for a given monomial. We are leaving this 



new tool for future work 1 1 
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Appendix A 



We present here an example of how the algorithm Simplif yAlllndex can be 
applied. The session is part of a Mathematica notebook. The inputs are indi- 
cated by In[] :=. The outputs have no prefix. The titles and comments are 
enclosed using the symbols (* comment*) . 

Since we use the Riemann tensor as an example we introduce this tensor 
using the TTC function InputSRiemann. This is internally equivalent to use 
InputTensor and InputSymmetries to declare the Riemann and Ricci tensors 
and to give the relation between them and between the curvature scalar. All 
tensor calculus conventions follow ( [fL2[ ). 

In this first short academic example we first introduce coordinate system 
cx4 and a metric gn together with the Riemann tensor. 

We define the tensor V and L with symmetries L l3k = L kjl , L y7c;m = L? lh < m 
and L^ k ' m + L lmj ' k + L lkm '^ = 0. Note that L have not separable properties so 
we will define all properties as a (2)-properties. 

Next we declare the indices for outputs and define and simplify a monomial 



In[] :=«ttc.m 



I Tools of Tensor Calculus 4.1.0 I 

I A.Balf agon.P.Castellvi and X.Jaen I 

I \protect\vrule widthOpt\protect\href {http : //balduf a. upc . es/ttc}{http : //balduf a.upc . es/ttc} I 

I e-mail :ttc@balduf a. upc.es I 

I version: September, 29,1999 I 



I Session started on I 

I October, 5 , 1999 I 

I at 11 h 15 min 48 s I 



In[] : = ( 

InputCoordinates [cx4,4] ; 
InputSMetric[gn,cx4,"g",g] ; 
InputSRiemann [gn, cx4, "R" , Rie ,Ric ,R] ; 
InputTensor [{L},cx4, {1,1,1}] ; 
InputTensor [{V},cx4,{l}] ; 
InputSymmetries [L [i , j ,k] , {i ,k} [2] ] ; 
InputSymmetries [L[i,j ,k, . ;m] , {i , j} [2] , 
Cyclic [j ,k,m] [2]] ; 
Input Index [{a,b, c ,d, e ,m,n}-] 
); 

3/2 L[i,j,k] Rie[-i,-j ,q,p] V[-k] //Index [gn] 

c d e a b 
3 L V R 
a b e c d 
+ 

2 

"/.//Index [gn , Simplif yAlllndex [2] ] 

b d c a 
3 L V R 
a b c d 
+ 

4 

b d c a 



13 



3 L V R 

c d 

+ 

4 

a deb 
3 L V R 

c d 



a d c b 

3 L V R 

c d 



a b 



3 L 



c d 



b a c d 

3 L V R 

c d 



4 

(* Note that if Ricci=0 then the original monomial 
is zero, which is not obvious from the beginning*) 

Next we show how the same kind of calculation is able to find a known result: 
The vanishing of the magnetic part of the Weyl tensor for static space-times, 
see ID 

(* first we introduce the predefined volume form which we name H*) 
In[]: = 

InputSAForm[gn,cx4, "H" ,H] ; 

(*Here we introduce the Killing field Vt, which we define as 
a vector antisymmetric between his own index and 
the covariant derivative index. *) 

In[] : = 

(InputTensor [Vt , cx4, {1}] ; 
InputSymmetries[Vt[a, . ;b] ,{{a,b}} [2]] ) ; 

(*Next we introduce the static condition. Since this 

property involves more than one tensor, it must be introduced 

explicitly by the user and taking into account 

its differential consequences. 

We introduce these conditions in BasicRules [2] 

using the specific syntaxis a_:>(a;b) instead of a_:>b !!*) 

(BasicRules [2] = Join [BasicRules [2] , { 

Vt[a_] Vt [b_ , . ; c_] : > (Vt [a] Vt[b,.;c];- 

Vt[c] Vt[a, . ;b]-Vt[b] Vt[c,.;a]), 

Vt[a_] Vt[b_, . ;c_, . ;d_] :>(Vt[a] Vt [b , . ; c , . ;d] ; 

-Vt[c] Vt[a, . ;b, . ;d]-Vt[b] Vt [c , . ; a, . ;d] - 

Vt[a,.;d] Vt[b, . ;c]-Vt[b, . ;d] Vt [c , . ; a] -Vt [c , . ; d] Vt[a,.;b]), 



Vt[b_,.;d_] Vt[c_, . ;a_] :>(Vt[b, . ;d] Vt[c,.;a]; 

-Vt[c] Vt[a, . ;b, . ;d]-Vt[b] Vt [c , . ; a, . ;d] - 

Vt[a] Vt[b, . ;c, . ;d]-Vt[a, . ;d] Vt [b, . ; c] -Vt [c , . ; d] Vt[a,.;b]), 
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}]); 

(* Next we define Weyl Tensor*) 
In[]: = 

IndexUpdate [Weyl ,":=", 

0[a,b,c,d]+Rie[a,b,c,d]-l/2 g[a,c] Ric[b,d] + 
1/2 g[a,d] Rlc[b,c]+l/2 g[b,c] Rlc[a,d]-l/2 g[b,d] Ric[a,c] + 
R/12(g[a,c] g[b,d]-g[b,c]g[a,d]+g[b,d]g[a,c]-g[a,d]g[b,c]) 
//Index [gn]] 

(*Next we define the magnetic part of Weyl, BW*) 
In[]: = 

IndexUpdate [BW," :=", 

1/2 H[a,r,c,m] Weyl [-c , -m,b,n] Vt [-r] Vt [-n] //Index [gn] ] ; 

(*We can see BW expanded. Now we like to use {i, j ,k,l. . .} for outputs*) 
In[] : =Input Index [{i , j , k, 1 ,m,n, o ,p}] ; 

In[] : =BW [-a, -b] //Index [gn, Super IndexExpand] 

lk lk lk 

R Vt Vt H R Vt Vt H Vt Vt H 

klij klij kli 

+ + 

i j 12 12 4 

1 m k 1 m k 1 

Vt Vt H R Vt Vt H R Vt Vt H 

kli mj klij m kli 

+ 

4 4 4 

1 m n k 
Vt Vt H R 

k 1 i m n j 

2 

(♦Finally we simplify BW, which give the final result*) 

In [] : ='/.//Index [gn, Simplif yAlllndex [2] ] 



(*Here we can see the internal aspect of some of the stored rules*) 
In [] : =L1 i st aMonomi s [gn , cx4 , 2 , 2 , 2] 

IndexD[l] IndexF[l] . ; IndexD [2] .; IndexD[3] 
{{Vt Vt 

IndexD [1] IndexD [2] IndexD [3] IndexF [2] 
H , {2, 1}}, 

IndexD [1] IndexD [2] IndexD [1] IndexD [3] IndexD [4] IndexF [1] 
{Vt Vt H 

IndexD [2] IndexF [2] IndexD [3] IndexD [4] 
R , {1, 2}}, 

In[] : = 

LlistaRules[gn,cx4,2,2,2] 

IndexD [1] IndexF [2] IndexD [2] IndexD [3] IndexD [4] IndexF [1] 
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{Vt Vt H 

IndexD[l] IndexD [2] IndexD[3] IndexD[4] 
R -> 0, 
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